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' Abstract. If neutrinos have masses and mixings, in general some will decay. In models 

^ |. with new particles and new couplings, some decay modes can be fast enough to be of 

' phenomenological interest. 

I INTRODUCTION 

(N 

>. 

I It is generally agreed that most probably neutrinos have non-zero masses and 

■ non-trivial mixings. This belief is based primarily on the evidence for neutrino 

mixings and oscillations from the data on atmospheric neutrinos and on solar neu- 



O i trinos [1]. 

O I If this is true, then in general the heavier neutrinos are expected to decay into 

the lighter ones via flavor changing processes. The only questions are (a) whether 



the lifetimes are short enough to be phenomenologically interesting (or are they 
, too long?) and (b) what are the dominant decay modes. 

Q I Throughout the following discussion, to be specific, I will assume that the neu- 

' trino masses are at most of order of a few eV. There are interesting things to be 

done with heavier decaying neutrinos which we will not discuss here today [2]. 

X 

S II RADIATIVE DECAYS 

For eV neutrinos, the only radiative decay modes possible are Ui —>■ uj + 7. They 
can occur at one loop level in SM (Standard Model). The decay rate is given by 
[3]: 



_9_ a Gl [Smtj 
16 vr 1287r3 



2.3 



rrii 



V f/* f/ 



ml 



m 



w , 



where 6mfj = mf — and a runs over e, /i and r. When rrii 3> mj, rrii ~ 0{eV) 
and for maximal mixing (4:11*^11^^) ~ 0(1) and {a = r) one obtains for F 



* Presented at the 5th International Conference on Physics Potential & Development of j-i^ ^ 
Colliders, San Francisco, December 1999. 



TsM ~ 10-^^ sec-^ 



(2) 



which is far too small to be interesting. If the dominant contribution is from 
electron (rather than r) then there is enhancement in matter (due to the electrons 
present) by a factor [4] 

(3) 



1024 (cc)-i J \m. 

where the function F tends to Arrii/Ei for relativistic u'-Si and this enhancement 
can be as large as 10^^. 

The decay mode i^i — > i^j + 'y comes from an effective coupling which can be 
written as: 

iPja^,{C + D-f5)il;i F^u (4) 



TTii + m 



Let us define kij as 



where /ib — e/2me. Since the experimental bounds on fj^^i, the magnetic moments 
of neutrinos, come from reactions such as i/ge — > e'V" which are not sensitive to 
the final state neutrinos; the bounds apply to both diagonal as well as transition 
magnetic moments and so can be used to limit and the corresponding lifetimes. 
The current bounds are [5]: 

/tg < 10-1° (6) 
k^ < 7.4.10-1° 
k^ < 5.4.10-^ 

For rrii ^ rrij, the decay rate for Ui i/j + ^ is given by 

This, in turn, gives indirect bounds on radiative decay lifetimes for 0{eV) neutrinos 
of: 

r^, > 5.10^^ sec (8) 
r^^ > 5.10^^ sec 
T^^ > 2.10^1 sec 



There is one caveat in deducing these bounds. Namely, the form factors C 
and D are evaluated at ~ 0{eV^) in the decay matrix elements whereas in the 
scattering from which the bounds are derived, they are evaluated at ~ 0{MeV'^). 
Thus, some extrapolation is necessary. It can be argued that, barring some bizarre 
behaviour, this is justified [6]. 

There are other bounds on radiative lifetimes of neutrinos which are all based 
on non-observation of the final state 7-ray. The most direct observational bounds 
from reactor and accelerators are [5]: 

T(i/e) > 300 s (9) 
t(i/^) > 15.4 s 

There are other more indirect bounds [5]. From cosmology: 

T > 2.10^^ s. (10) 

Prom x-ray and 7-ray fluxes: 

T > T.IO'' s. (11) 

Here we have set rrii ~ 0{eV). All these bounds depend on assuming that rrii ^ rrij 
in the mode Ui — > i/j + 7. The bounds are not valid if there is near degeneracy; 
since in this case the 7-ray can be soft. 



Ill THE SCIAMA MODEL 

Dennis Sciama proposed an intriguing and imaginative scenario for radiative 
decay of a neutrino about 10 years ago [7]. It went thru several metamorphoses; 
starting out as a decay mode i't ^ + 1 ^-nd ending up most recently as i/gt 
+ 1 [8]. 

The proposal was this: (a) suppose the decay i^a ^ + 1 takes place, (b) the 
mass of Uairriiya) is ~ 27.4 eV, and m^^ ^ m,^^, (c) thus 13.7 eV for a slow 

moving (d) the lifetime by this mode for i/^ is about 2.10^^ sec. These properties 
are consistent with all the bounds. The purpose was to account for the anomalous 
amount of ionised H in interstellar space. The energy of the decay photon is just 
right to ionise hyrogen; and the lifetime is chosen so as to provide the required 
amount of ionization. 

After a lengthy innings (of 10 years), it appears that the Sciama model is finally 
ruled out. The most recent evidence against it is the non-observation of the line at 
911 A (corresponding to 13.7 eV) with the predicted intensity [9]. 

IV INVISIBLE DECAYS 

A decay mode with essentially invisible final states which does not involve any 
new particles is the three body neutrino decay mode. The decay is Ui —>■ VjUjUj. 
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This decay, like the radiative mode, can occur at one loop level in SM. With a mass 
pattern mj ^ nij the decay rate can be written as 

^ (12) 



1927r3 

In the SM at one loop level, with the internal r dominating, the value of is given 
by [10] 



4m-^ Z ^ (Ur,U-f (13) 




With maximal mixing ~ 3.10" . Even if e were as large as 1 with new 
physics contributions; it only gives a value for F of 5.10~^^ sec~^. Hence, this 
decay mode will not yield decay rates large enough to be of interest. It is worth 
noting that the current experimental bound on e is quite poor: e < O(IOO) [11]. 



If the only new particle introduced is an 1=0, L=0, J=0, massless particle such 
as a Goldstone boson (familon), then a new decay mode is possible: 

^c^L ^^Pl + X (14) 
which arises from a coupling of the form 

Qp^pLlu.'^aLdi.X (15) 

with a decay rate 

r = ^ (16) 

SU(2)i symmetry predicts a similar coupling for the charged leptons with decay 
mode — ^ + X- Thus the Va lifetime is related to the B.R. (£„ (-p + x) by 



B.R.{(i^-^f.p + X) \mej 
The current bounds on /i and r branching ratios are [5,12] 

B.R.in^ex) <2.10-^ (18) 
B.R.{t^HX) <7.10-^ 

and lead to 

r., > 10^^ s (19) 
r.. > 10^° s. 

Coupling to a scalar field with 7 = 1/2, L = would have the same constraint 
as above; and in addition would require fine tuning to avoid mixing with the SM 
Higgs. 
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VI FAST INVISIBLE DECAYS 



The only possibility for fast invisible decays of neutrinos seems to lie with Ma- 
joron models. There are two classes of models; the 1=1 Gelmini-Roncadelli [13] 
majoron and the 1=0 Chikasige-Mohapatra-Peccei [14] majoron. In general, one 
can choose the majoron to be a mixture of the two; furthermore the coupling can 
be to flavor as well as sterile neutrinos. The effective interaction is of the form: 

gai^pi'a J (20) 

giving rise to decay: 

Ua^J^p +J (21) 

where J is a massless J = L = 2 particle; and are mass eigenstates which 
may be mixtures of flavor and sterile neutrinos. Models of this kind which can 
give rise to fast neutrino decays and satisfy the bounds below have been discussed 
by Valle, Joshipura and others [15]. These models are unconstrained by fi and r 
decays which do not arise due to the AL = 2 nature of the coupling. The 1=1 
coupling is constrained by the bound on the invisible Z width [2]; and requires that 
the Majoron be a mixture of 1=1 and 1=0. The couplings of ly^ and {g^ and Qe) 
are constrained by the limits on multi-body tt, k decays vr — > jii/i/v and K — > jii/vi/ 
and on /i — e university violation in tt and K decays [16]. 

Granting that models with fast, invisible decays of neutrinos can be constructed, 
can such decay modes be responsible for any observed neutrino anomaly? 

We assume a component of i.e., U2, to be the only unstable state, with a 
rest-frame lifetime Tq, and we assume two flavor mixing, for simplicity: 

Va — cos6v2 -\- sinOvi (22) 

with 777.2 > rni- From Eq. (2) with an unstable 1^2, the Va survival probability is 

Paa = sin^d + cos^eex.^{-aL/E) (23) 
+ 2sin'^ecos^eex-p{-aL/2E)cos{5m^L/2E), 

where Sm^ — ml — ml and a = 7772/To. Since we are attempting to explain neutrino 
data without oscillations there are two appropriate limits of interest. One is when 
the Sm"^ is so large that the cosine term averages to 0. Then the survival probability 
becomes 

P^^^sin^e +cos'^0exp{-aL/E) (24) 

Let this be called decay scenario A. The other possibility is when ^777^ is so small 
that the cosine term is 1, leading to a survival probability of 
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(25) 



corresponding to decay scenario B. Decay models for both kinds of scenarios can 
be constructed; although they require fine tuning and are not particularly elegant. 

The possibility of solar neutrinos decaying to explain the discrepancy is a very 
old suggestion [17]. The most recent analysis of the current solar neutrino data 
finds that no good fit can be found: Uei ~ 0.6 and Ti, (E=10 MeV) ~ 6 to 27 
sec. come closest [18]. The fits become acceptable only if the suppression of the 
solar neutrinos is energy independent as proposed by several authors [19] (which is 
possible if the Homestake data are excluded from the fit). The above conclusions 
are valid for both the decay scenarios A as well as B. 

For atmospheric neutrinos, assuming neutrino decay scenario A, it was found 
that it is possible to choose 9 and a to provide a good fit to the Super-Kamiokande 
L/E distributions of events and z/^/^e event ratio [20]. The best-fit values of the 
two parameters are cos'^9 0.87 and a ~ IGeV/ D^. where De = 12800 km is the 
diameter of the Earth. This best-fit a value corresponds to a rest-frame 1^2 lifetime 
of 

To - m2/a ~ X 10-i°s. (26) 

(ley) 

However, it was then shown that the fit to the higher energy events in Super-K 
(especially the upcoming muons) is quite poor [21]. The reason that the inclusion 
of high energy upcoming muon events makes the fits poorer is very simple. In the 
decay A scenario, due to time dilation, the decay is suppressed at high energy and 
there is hardly any depletion of v'^s, contrary to observation. 

Turning to decay scenario B, consider the following possibility [22]. The three 
states i^^, i^r, J^s (where Us is a sterile neutrino) are related to the mass eigenstates 
1/21^31^4 by the approximate mixing matrix. 



(27) 



and the decay is 1^2 ^ 1^4 + J. The electron neutrino, which we identify with ui, 
cannot mix very much with the other three because of the more stringent bounds 
on its couplings [16], and thus our preferred solution for solar neutrinos would be 
small angle matter oscillations. 

In this case the Sml^ in Eq. (23) is not related to the 5m|4 in the decay, and 
can be very small, say < 10~^eV^ (to ensure that oscillations play no role in the 
atmospheric neutrinos). Then the oscillating term is 1 and Piiy^ — > f^) is given by 
Eq. (25). 

In order to compare the predictions of this model with the standard Ur 
oscillation model, we have calculated with Monte Carlo methods the event rates 
for contained, semi-contained and upward-going (passing and stopping) muons in 

the Super-K detector, in the absence of 'new physics', and modifying the muon 
neutrino fiux according to the decay or oscillation probabilities discussed above. 
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We have then compared our predictions with the SuperK data [23], calculating 
a to quantify the agreement (or disagreement) between data and calculations. 

In performing our fit (see Ref. [22] for details) we do not take into account any 
systematic uncertainty, but we allow the absolute flux normalization to vary as a 
free parameter /5. 

The decay model of Equation (24) above gives a very good fit with a minimum 
— 33.7 (32 d.o.f.) for the choice of parameters 



and normalization (3 = 1.17. 

In Fig. 1 we compare the best fits of the two models considered (oscillations and 
decay) with the SuperK data. In the figure we show (as data points with statistical 
error bars) the ratios between the SuperK data and the Monte Carlo predictions 
calculated in the absence of oscillations or other form of 'new physics' beyond the 
standard model. In the six panels we show separately the data on e-like and fi- 
like events in the sub-GeV and multi-GeV samples, and on stopping and passing 
upward-going muon events. The solid (dashed) histograms correspond to the best 
fits for the decay model <-> u-r oscillations). One can see that the best fits of 
the two models are of comparable quality. The reason for the similarity of the 
results obtained in the two models can be understood by looking at Fig. 2, where 
we show the survival probability P{i'n I'^i) of muon neutrinos as a function of 
L/E^ for the two models using the best fit parameters. In the case of the neutrino 
decay model (thick curve) the probability P{i^^ i^^) monotonically decreases 
from unity to an asymptotic value sin*^^ ~ 0.49. In the case of oscillations the 
probability has a sinusoidal behaviour in L/E^. The two functional forms seem 
very different; however, taking into account the resolution in L/E^, the two forms 
are hardly distinguishable. In fact, in the large L/E^, region, the oscillations are 
averaged out and the survival probability there can be well approximated with 0.5 
(for maximal mixing). In the region of small L/Ey both probabilities approach 
unity. In the region L/Ey around 400 km/GeV, where the probability for the 
neutrino oscillation model has the first minimum, the two curves are most easily 
distinguishable, at least in principle. 

Decay Model 

There are two decay possibilities that can be considered: (a) V2 decays to z/4 
which is dominantly Vg with V2 and mixtures of and z/^, as in Eq. (27), and 
(b) 1^2 decays into which is dominantly z/^ and V2 and are mixtures of and 
Vg. In both cases the decay interaction has to be of the form 



where J is a Majoron field that is dominantly iso-singlet (this avoids any confiict 

with the invisible width of the Z). Viable models for both the above cases can 
be constructed [15,24]. However, case (b) needs additional iso-triplet light scalars 
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Ty/rriy = 63 km/GeV, cos^ 9 = 0.30 



(28) 



(29) 



which cause potential problems with Big Bang Nucleosynthesis (BBN), and there 
is some preliminary evidence from Super K against Vfx-i's mixing [25]. Hence we 
only consider case (a), i.e. z/2 — * z/4 + J with Vs, as implicit in Eq. (27). With 
this interaction, the U2 rest-lifetime is given by 

where 5m^ = m| — ml and x = 777.4/777.2 {0 < x < 1). Prom the value of — 
T2/7772 = 63 km/GeV found in the fit and for x = 0, we have 

/5m2~0.16eV^ (31) 

Combining this with the bound on g"^ from K ^ /j, decays of g^ < 2.4 x 10~"^ [16] 
we have 

Sm^ > 650 eV^ . (32) 

Even with a generous interpretation of the uncertainties in the fit, this SmF' implies 
a minimum mass difference in the range of about 25 eV. Then 1/2 and are nearly 
degenerate with masses > (9(25 eV) and z/4 is relatively light. We assume that a 
similar coupling of to 2/4 and J is somewhat weaker leading to a significantly 
longer lifetime for 1/3, and the instability of is irrelevant for the analysis of the 
atmospheric neutrino data. 

For the atmospheric neutrinos in SuperK, two kinds of tests have been proposed 
to distinguish between t'^-^'r oscillations and u^-Ug oscillations. One is based on the 
fact that matter effects are present for p^-Ug oscillations [26] but are nearly absent 
for lyn-Vr oscillations [27] leading to differences in the zenith angle distributions due 
to matter effects on upgoing neutrinos [28]. The other is the fact that the neutral 
current rate will be affected in z/^-z/g oscillations but not for v^-Vr oscillations as 
can be measured in events with single 7r°'s [29]. In these tests our decay scenario 
will behave as a hybrid in that there is no matter effect but there is some effect in 
neutral current rates. 

Long-Baseline Experiments 

The survival probability of function of L/E is given in Eq. (1). The 

conversion probability into Ur is given by 

P{^^, ^ ^r) = sin^ ^cos^ ^(1 - e-'^'^l'^^f . (33) 

This result differs from 1—P{v^ — > i/^) and hence is different from v^-Vt oscillations. 
Furthermore, P^i^n y^) + P{i'i^ — > 1/-,-) is not 1 but is given by 

P{i^i. ^ i^m) + Pi^^^ ^iyr) = l- cos^ 9(1 - e-"^/^) (34) 

and determines the amount by which the predicted neutral-current rates are affected 
compared to the no oscillations (or the t'^-i^r oscillations) case. In Fig. 3 we give 



the results for Piv/j, v^), -P(^'j^ i't) and P{yfx — > v^) + — ^ i^t) for 
the decay model and compare them to the Vh-Vt oscillations, for both the K2K 
[30] and MINOS [31] (or the corresponding European project [32]) long-baseline 
experiments, with the oscillation and decay parameters as determined in the fits 
above. 

The K2K experiment, already underway, has a low energy beam Ejy fa 1-2 GeV 
and a baseline L = 250 km. The MINOS experiment will have 3 different beams, 
with average energies E,^ = 3, 6 and 12 GeV and a baseline L = 732 km. The 
approximate L/ E^, ranges are thus 125-250 km/GeV for K2K and 50-250 km/GeV 
for MINOS. The comparisons in Figure 3 show that the energy dependence of 
survival probability and the neutral current rate can both distinguish between the 
decay and the oscillation models. MINOS and the European project may also have 
r detection capabilities that would allow additional tests. 

Big Bang Nucleosynthesis 

The decay of V2 is sufficiently fast that all the neutrinos (z/g, ^'^t, '^t, t-'s) and the 
Majoron may be expected to equilibrate in the early universe before the primordial 
neutrinos decouple. When they achieve thermal equilibrium each Major ana neu- 
trino contributes = 1 and the Majoron contributes = 4/7 [33], giving and 
effective number of light neutrinos N^, = 4| at the time of Big Bang Nucleosynthe- 
sis. From the observed primordial abundances of "^He and ^Li, upper limits on N^, 
are inferred, but these depend on which data are used [34-36]. Conservatively, the 
upper limit to N^, could extend up to 5.3 (or even to 6 if ^Li is depleted in halo 
stars [34]). 

Cosmic Neutrino Fluxes 

Since we expect both U2 and z/3 to decay, neutrino beams from distant sources 
(such as Supernovae, active galactic nuclei and gamma-ray bursters) should contain 
only i/g and but no Ufj,, u^, Vt and Vr- This is a very strong prediction of our 
decay scenario. We can compare the very different expectations for neutrino flavor 
mixes from very distant sources such as AGN's or GRB's. Let us suppose that at 
the source the flux ratios are typical of a beam dump, a reasonable assumption: 
A^zA. • '^Vf, '■ = 1:2:0. Then, for the conventional oscillation scenario, when all 
the 5m^'fi satisfy 5rn? L/AE » 1, it turns out curiously enough that for a wide 
variety of choices of neutrino mixing matrices, the final fiavor mix is the same, 
namely: N^,^ : N^,^ : Ny^ = 1:1:1. In the case of the decay B scenario, as 
mentioned here, we have N^^ : Ny^ : Ny^ = 1:0:0. The two are quite distinct. 
Techniques for determining these flavor mixes in future KM3 neutrino telescopes 
have been proposed [37]. 

Reactor and Accelerator Limits 

The Ve is essentially decoupled from the decay state V2 so the null observations 
from the CHOOZ reactor are satisfied [38] . The mixings of t'^ and Vt with fg and 
Uf. are very small, so there is no confiict with stringent accelerator limits on fiavor 
oscillations with large 5w? [39]. 
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Summary 

In summary, neutrino decay remains a viable alternative to neutrino oscillations 
as an explanation of the atmospheric neutrino anomaly. The model consists of two 
nearly degenerate mass eigenstates z/2, z/3 with mass separation > (9(25 eV) from 
another nearly degenerate pair i/i, 1/4. The and i^r flavors are approximately 
composed of 1/2 and 1/3, with a mixing angle 623 — 57°. The state 1^2 is unstable, 
decaying to P4 and a Majoron with a hfetime T2 ~ 10~^^ sec. The electron neutrino 
z/g and a sterile neutrino Ug have negligible mixing with z/^, z/^ and are approximate 
mass eigenstates (z/g ~ ^i-^s ~ ^^4)1 with a small mixing angle 614 and a 6ml^ ^ 
10"^ eV^ to explain the solar neutrino anomaly. The states z/3 and z/4 are also 
unstable, but with 1/3 lifetime somewhat longer and 1/4 lifetime much longer than the 
z/2 lifetime. This decay scenario is difficult to distinguish from oscillations because 
of the smearing in both L and E„ in atmospheric neutrino events. However, long- 
baseline experiments, where L is fixed, should be able to establish whether the 
dependence oi L/ is exponential or sinusoidal. In this scenario only ui is stable. 
Thus, neutrinos of supernovae or of extra galactic origin would be almost entirely 
z/g. The contribution of the electron neutrinos and the Majorons to the cosmo logical 
energy density VL is negligible and not relevant for large scale structure formation. 

VII CONCLUSION 

In general, when neutrinos have masses and do mix; they can decay as well as 
oscillate. Some neutrino anomalies may be caused by one or the other or both. 
Non-oscillatory explanations have to be ruled out experimentally. Eventually, in- 
visible neutrino decays can be severely constrained when data from Long Baseline 
experiments with good resolution in L and E become available. 
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FIGURE 1. Comparison of decay model (solid histograms) and v^-i't oscillation model (dashed 
histograms) with SuperK data. 
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FIGURE 2. Survival probabiliity for j/^ versus \ogiQ{L/E) for the decay model (heavy solid 
curve) and oscillation model (thin curve). 
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FIGURE 3. Long-baseline expectations for the K2K and MINOS long-baseline experiments 
from the decay model and the v^-Vr oscillation model. The upper panel gives the neutral current 
predictions compared to no oscillations (or v^-Vr oscillations). 
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